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I Abstract. We discuss and construct an action functional (logarithmic action) for the 

simply connected Laplacian growth and obtain its variation. This variation admits various 
interpretations. In particular, we consider a general smooth subordination evolution and 
give connections with the Virasoro algebra and Neretin polynomials. 
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*0 ; 1. Introduction 
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An important quantity in mechanics and the field theory is the action, a local functional 
over maps <fi between the space-time £ and the target space which we assume to be the 
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Euclidean straight line M.. For example, one may consider the action given by the Dirichlet 
integral 

S[<j>] = J \V<J>\ 2 dv, 



E 



where the metric structure of S and the volume element dv are to be taken into account as 
well as certain properties of smoothness of <j), that give sense to the right-hand side of the 
above equality. In the classical setup of mechanics this action represents the energy of the 
system as infinitesimally it is just the scalar product of the field V0 and its momentum 
V(f)dv for the potential 0. We note that in general, the action is the integral over the 
classical Lagrangian. If time does not enter explicitly into the Lagrangian, then the system 
is closed, and a typical example of such a Lagrangian is the kinetic energy minus the po- 
tential energy. One may assume different smooth functions on TS as Lagrangians, however 
physical or geometrical background of the underlying space £ forces certain restrictions in 
the choice. Nevertheless, different points of view on the same object can lead to different 
functionals as Lagrangians. 

The classical field theory studies the extremum of the action functional, and its critical 
value is called the classical action. The critical point <p* satisfies Hamilton's principle (or 
the principle of the least action), i.e., &S[0*] = 0, which is the Euler-Lagrange equation 
for the variational problem defined by the action functional. For the action given by the 



2000 Mathematics Subject Classification. Primary: 76D27, 81T40; Secondary: 30C35, 30C62, 30F45, 
30F60. 

Key words and phrases. Liouville action, Laplacian growth, Hele-Shaw problem, conformal map, Vira- 
soro algebra. 

1 



2 



ALEXANDER VASIL'EV 



Dirichlet integral the classical action is achieved for the harmonic 0* and the principle of 
the least action leads to the Laplacian equation A0 = 0. 

The Liouville action plays a key role in two-dimensional gravity. It is based on a La- 
grangian given on the base of the Riemannian geometry of the underlying space. Liouville 
action describes a highly non-trivial dynamics in quantum field theory and appears in con- 
nection with Feynman's path integral that represents the transition amplitude between two 
quantum states of a system expressed as a sum over contributions from possible classical 
histories of that system. Complex transition functions appear naturally in the theory of 
evolution equations. The original formulation of quantum Liouville theory through path 
integral has been obtained by Polyakov (SHj in 1981 where the domain of integration con- 
sisted of all smooth conformal metrics ds 2 on an n-punctured Riemann sphere (n > 4 to 
guarantee the hyperbolicity) . A thorough mathematical treatment has been made later by 
Takhtajan and Zograf (see (HE], US], jH, H3). 

It is not very surprising that several "quantum features" appear in non-linear problems 
of hydrodynamics, in particular, in the Laplacian Growth problem. In 1898 Hele-Shaw (THj 
proposed his famous cell that was a device for investigating a flow of viscous fluid in a 
narrow gap between two parallel plates. 

The dimensionless model of a moving viscous incompressible fluid in the Hele-Shaw cell 
is described by a potential flow with the velocity field V = (Vi, V 2 ). The pressure p is the 
potential for the fluid velocity 

h 2 

v = ~ik Vp ' 

where h is the cell gap and \i is the viscosity of the fluid (see, e.g. [27J Ell)- 

Given an incompressible advancing fluid injected through a point source, the Laplacian 
growth is formulated as a moving boundary problem for the Laplacian equation for a 
function p(z,t), supported in a domain Q(t) C C as a function of z, 

(1) Ap = -2tt5 (z), z e n(t), 

where t is the time variable and do is the Dirac distribution supported in 0. The dynamic 
boundary condition is given by putting 



(2) p 



= 0, 

8Q(t) 



and the kinematic condition for the motion of the boundary dQ(t) is given by the normal 
velocity as 

dp 

(3) v n = -—, x e dn(t), 

where n is the outward unit vector to dQ(t). 
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Through the similarity in the governing equations, Hele-Shaw flows can be used to study 
models of saturated flows in porous media governed by Darcy's law. Over the years various 
particular cases of such a flow have been considered. Different driving mechanisms were 
employed, such as surface tension or external forces (suction, injection). We mention here 
a 600-paper bibliography of free and moving boundary problems for Hele-Shaw and Stokes 
flows since 1898 up to 1998 collected by Gillow and Howison jH]. 

As it has been shown in 1221 123 HI] , the Laplacian growth problem can be embed- 
ded into a larger hierarchy of domain variations (Whitham-Toda hierarchy) for which all 
Richardson's complex moments [33] are treated as independent variables (generalized time 
variables), and form an integrable system. Finally, the Laplacian growth has been modeled 
in the moduli space of Riemann surfaces |21j . 

In the classical simply connected case of the Laplacian growth without gravity the kinetic 
energy is given by the Dirichlet integral for the pressure as a potential. However, given 
the evolution of the phase domain in time as a closed system with a Riemannian metric as 
the geometrical background, one may construct the action functional based on a different 
Lagrangian. The idea of the construction of the Liouville action gives us a way to derive 
an action functional for the Laplacian growth. The aim of this paper is to construct the 
logarithmic action for the Laplacian growth and to obtain its variation. Then we shall 
study a general smooth subordination and interpret the variation of the logarithmic action 
through the infinitesimal version of the action of the Virasoro-Bott group on the space of 
analytic univalent functions. 

The author thanks Professor Leon Takhtajan for discussions on this subject. 



2. Liouville and logarithmic actions 

Following the classical approach by Poicare [2H] we consider a compact Riemann surface 
S of genus g > 2 that admits uniformization by a Fuchsian group G acting in the unit disk 
U, S = U/G. If z is a local analytic coordinate defined on an open set, then the Riemann 
metric ds 2 is represented as ds 2 = p 2 \dz\ 2 for a density p. A conformal metric corresponds 
to the invariancy under the change of the local coordinate. If z' is another local coordinate 
defined in an open set and ds 2 = (p') 2 \dz'\ 2 in terms of the local coordinate z', then we 
have p(z) = p' (z' (z))\dz' / dz\ in the intersection of these sets. 

The Gaussian (sectional) curvature x of this metric is calculated by the formula x = 
— A log p. It follows from the uniformization theorem that there exists a unique conformal 
metric of Gaussian curvature x = — 1 which is called the Poicare (hyperbolic) metric. In 
terms of the parameter z this is equivalent to the Liouville equation 
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where p = log p 2 . In the simplest case of the unit disk the Poicare metric is given as 

4KI 2 



ds z 



:i-ici 



2\2 ' 



Considering the universal covering of S by the unit disk and the automorphic (with respect 
to G) projection h : U — > S we deduce that 

ds 2 



A\h- V 


) 2 \d(\ 2 




h- 1 ] 


2\2 



where ( is taken from the fundamental polygon U/G. This metric is complete and the area 
of S in this metric is 47r(g — 1) (by the Gauss-Bonnet theorem). 

It seems that equation (0J is the Euler-Lagrange equation for the variational problem 
defined by the functional 

(5) S[<j>] = j (|0,| 2 + e+)da z , 

U/G 

where da z = I dzh 2 dz which may be chosen as the Liouville action. However, the deal is 
much more difficult, because the integrant (the first term) does not keep the local 2- form 
\4> z \ 2 dz A dz invariant under the change of the local parameter from chart to chart on 
the Riemann surface. Thus, the functional (j5J) is well-defined only for simply connected 
domains and the correct Liouville action requires an additional term which has been given 
by Takhtajan and Teo in |37j . 

When the underlying Riemann surface have singularities (e.g., punctures, branch points) 
the action functional is not well defined globally on the surface either. In this case, there 
are two options to treat the problem. One of them is to change the metric, the other is 
to make certain regularization. Actually, Poincare [2H] worked in the presence of parabolic 
singularities. In this case some additional terms appear in (|SJ), see [13]. But certainly one 
may look for equilibrium between these two options. 

Let us turn now to the Laplacian growth. As it was said in Introduction, the problem 
(dHHl) defines the Laplacian growth. More thoroughly we give a strong formulation of this 
problem. Let Q(t) cC, 6 &{t), be a one-parameter family of bounded domains. We call 
the family {fi(£)} smooth if dQ(t) are smooth (C°°) interfaces for each t, and the normal 
velocity v n continuously depends on t at any point of dfl(t). Each Q(t) is supposed to be 
simply connected for any t G [0,T) fixed. A smooth family Q(t), < t < T, as above, 
is said to be a strong solution for the Laplacian growth if there exists a potential p(z, t), 
z G fl(t), such that all conditions (PQ-EJ) a re satisfied. The family of Q(t) forms a strong 
subordination chain of bounded domains: Q(s) C Q(t) for < s < t < T, 6 Q(0). 
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It is known that if the initial domain Q(0) has an analytic smooth boundary, then the 
strong solution exists locally in time until the boundary dQ(t) develops a cusp in a blow-up 
time or Q(t) changes its topology (see, e.g. [TBJ E21 GH1 EE])- 

By the Riemann theorem we construct a conformal time- dependent map z = f{(,t) 
from the unit disk U onto the phase domain Q(t), f(0,t) = 0, f'(0,t) > 0. The function 
/(C) 0) = /o(C) parameterizes the initial boundary dQ(0) = {fo(e w ),9 e [0, 2n)} and the 
moving boundary is parameterized by dQ(t) = {f(e te , t),9 e [0, We use the notations 

f = df/dt,f' = df/d(. 

Let us consider the complex potential W(z,t), Re W = p. For each fixed t it is a 
multivalued analytic function defined in Q(t) whose real part solves the Dirichlet problem 
Making use of the Cauchy- Riemann conditions we deduce that 

dW dp dp 

= a *^~> z = x + iy. 

a z ox ay 

Since Green's function solves (JTJ^, we have the representation 
(6) W(z,t) = -log z + w (z,t), 

where wo(z,t) is an analytic regular function in Q(t). Because of the conformal invariance 
of Green's function we have the superposition 

(Wo /)(£*) = -log C, 



and the conformally invariant complex velocity is just W'(z,t) = —jrrr(z,t), where £ = 
f~ l {z,t) is the inverse to our parametric function / and prime means the complex deriva- 
tive. Rewriting this relation we get 



(7) (W'(z,t)dzf = ^. 



The velocity field (— Vp) is the conjugation of (— W). In other words the velocity field is 
directed along the trajectories of the quadratic differential in the left-hand side of (J7J) for 
each fixed moment t. The equality implies that the boundary dQ(t) is the orthogonal 
trajectory of the differential (W'(z, t) dz) 2 with a double pole at the origin. The dependence 
on t yields that the trajectory structure of this differential changes in time, and in general, 
the stream lines are not inherited in time. These lines are geodesic in the conformal metric 
\W'(z, t)\\dz\ generated by this differential. 

Observe that being thought of as a Riemann surface, the phase domain Q(t) \ {0} is 
hyperbolic and it admits the Poincare metric with constant negative curvature 



If- 1 '] 


2 


I/- 1 ! 


2 log 2 If' 1 ] 



The asymptotics about the origin and close to the hyperbolic boundary implies that the 
standard expression (jSJ) for the Liouville action can not be used any more. Moreover, 
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the action integral for the hyperbolic metric in the punctured unit disk has the following 
asymptotics 



J (\<p z \ 2 + e*)da z ~ -27rlog£i - 4vr log | log e x | 



4tt 

+ 47rlog| log(l - e 2 ) 



log(l - e 2 ) 

where U e — {z : E\ < \z\ < 1 — 82}- Therefore, the corresponding integral for Q(t) \ {0} 
will have a similar asymptotics plus terms containing f'(0,t) and the boundary distortion 
(in our case the boundary derivatives) at the unit circle that makes it difficult to operate 
with. 

Let us use the flat logarithmic metric instead generated by ((7j) which seems to be more 
natural for the Laplacian growth 

d^ = ^^\dz\ a =\W f \ 2 \dz\ 3 . 

The hyperbolic boundary is not singular for this metric whereas the origin is. But it is a 
parabolic singularity which can be easily regularized. 

The density of this metric satisfies the usual Laplacian equation ip zz = in Q(t) \ {0}, 

I f — 1 ' I 2 

where (p(z) = log jj=rp ■ Obviously the Laplacian equation is the Euler- Lagrange equation 
for the variational problem defined by the Dirichlet integral 

2 da z , 



D 

locally for any measurable set D C Q(t) \ {0}. However, this functional cannot be defined 
globally in Q(t) \ {0} because of the parabolic singularity at the origin. To overcome this 
obstacle we define the classical action in the following way. Let fl £ (t) = Q(t) \ {z : \z\ < e} 
for a sufficiently small s, U e = {( : e < |£| < 1}. The function ip possesses the asymptotics 

ip ~ log \^ z \ ~ — - as z -> 0. 



1 19 ' 


\<fz\ ~ TT 


\ Z \ 


\z\ 



Therefore, the finite limit 

(8) limi / \ip z \ 2 da z + 2Tr\ogs } =: S[<p] 




exists and we call it the logarithmic action for the Laplacian growth. 



Lemma 1. The Euler-Lagrange equation for the variational problem for the logarithmic 
action S[<p} is the Laplacian equation A0 = — 47r<5o(^), z G Q(t), where Sq(z) is the Dirac 
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distribution supported at the origin, where <fi is taken from the class of twice differentiable 
functions in Q(t) \ {0} with the asymptotics <p ~ — log \z\ 2 as z — > 0. 

Proof. Let us consider first the integral 



J Wz\ 2 da z = J Xn £ (t)\fz\ 2 da z 



n E (t) 



where Xn s (t) is the characteristic function of £l e (t). Then, due to Green's theorem, 



lim 



S E [<j> + hu]-S E [4>] 
h 



2 J KM** **>.*>. 



(9) 



uA(f)da z -\ — / u—ds 



dn 



n e (t) 



an e (t) 



in distributional sense for every C°°(C) test function u supported in fl(t). On the other 
hand, we have d<p/dn ~ — 2/e as e — > and « = on 90(t). Therefore, the expression © 
tends to 

uA(j)do~ z — 2nu(0), 



n(t) 

as e — > 0, and the latter must vanish, that is equivalent to the Laplacian equation mentioned 
in the statement of the lemma. Obviously, the logarithmic term in the definition of S[<p] 
does not contribute into the variation. □ 

Straightforward calculation gives 
Hence, the action S can be expressed in terms of the parametric function as 

2 



1 ff " 7 +-,)or\z,t.) 



(10) 



S[p] = S[f] = lim J 



rfa c + 27rlog £ V+27rlog|/ , (0,t)|, 



or adding the logarithmic term into the integral we obtain 
ID S\f\ 



u 



r i 

/' c 



- |4)^ C + 27Tlog|/ / (0,t)|. 



The functional (|TT|) resembles the universal Liouville action defined by Takhtajan and Teo 
in [3B1 IB9J for quasicircles which is based on conformal maps from the unit disk and from 
its exterior. 
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Observe that the classical kinetic energy £ for the harmonic potential p is calculated by 
the Dirichlet integral 



\Pz\ 2 da z = / \W'\ 2 da 



D D 

locally for any measurable set D C Q(t) \ {0}. However, this integral again cannot be 
defined globally in Q(t) \ {0}. Treating £ in the same way as S we come to the following 
finite limit 

£:=lim< J \W'\ 2 d(j z + 2ii\oge 

or in terms of the parametric function / 

£ = £[/] = 27rlog|/'(0,f)|. 

The latter expression allows us to think of £ as a capacity which exactly corresponds to 
the physical sense of £ as minimal energy. 

3. Variation of the logarithmic action 

The Laplacian growth problem (JIHsIJ) being rewritten for the parametric time dependent 
function /:[/—> Q(t) admits the form of the so-called Polubarinova-Galina equation, which 
is in principle, a reformulation of the kinematic condition (see, e.g., |16[ I46j). Namely, the 
function f((,t) satisfies the non-linear first-order partial differential equation 



(12) Re [f Cf) = 1, ICI = 1, 

with the initial condition /(CO) = /o(C)- We denote by Sf the Schwarzian derivative 

f" 3ff"\ 2 

Sf -y-2{j) ' 

and by 

Re (l + ^ 



\f'(*°,t)\ 

the curvature of the boundary dQ(t) at the point f(e l9 ,t). 

Theorem 1. Let z = f{(,t) be the parametric function for the Laplacian growth, £[f] be 
the kinetic energy, and S[f] be the logarithmic action. Then 

2tt 
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and 



d 

dt 



2tt 



2tt 



(S[f\ + £ If]) = 2 / x\6, t)M + 



MPW* 8 {e2 " s,)de - 



Proof. Making use of the Cauchy-Schwarz representation we extend this equation into the 
unit disk 



(13) 
where 

(14) 



f = (f'p((,t), 



2tt 



p(C,t) 



1 e w + ( 
2tt J \f'{e i9 ,t)\ 2 e ie -( 



d9. 



Immediately, we obtain that 



|£[/] = |2.1og|/(0, f )| 



2tt 



|/'(e",t)|-- 



:d6. 



and hence, 



2?r 

(15) ^[/]=2Re |(T + ^ (1 + C ^ Wf , ()+ ^ (f , () )^ + |_l_ (i9 . 



The integral in the first term of the right-hand side of (fl3j) we rewrite by Green's theorem 
as 



2i 



re 



i + Cj)p(C,t) + Cp'(Ct)) d(- 



2- 



2 J |/'(e*V)| 2 



de, s 1 = au, 



taking into account a singularity at the origin. Applying the Cauchy-Schwarz formula to 
the first term in I containing p we arrive at 



2tt 



2 Re J 



l + e* 



/'(e<«t) 



da 



|/'(e* t)| 2 



2- 



2t 



+ Re / 1 + e 



f'(e ia ,t) 



e i y(e tQ ,t)da- 



\f'(e id ,t)\ 



-d9, 
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or 



2tt 



2tt 



t)\ 



-de 



■ f" 

1 + e ta — 

r 



da 



2vr 



+ 



2tt 

/ im (i^g^w)^. 



These equalities are thought of as limiting values making use of the analyticity of / on the 
boundary and f'((,t) ^ for all ( in the closure of the unit disk. Let us denote in the 
latter expression by J 2 the last integral and by J± the intermediate one. We have, 



J 2 



2vr 

/ Im ( 



1 + e iaJ — I Im 



2-k 



2e'V 



Obviously, 

d fe ie + C\ ~2(ie l 



f'J 2ir J \f'(e ie : t)\ 2 (e ie -e 







ia\2 



d9 da. 



-, and — 



2 in*"?" 



d0\e ie -(J (e^-C) 2 ' 09 \\f'(e ie ,t)\ 2 J \f'(e id ,t)\ 2 f ' 

Integrating by parts and applying the Cauchy-Schwarz formula again we obtain 

2tt 



o 



dd 



f'JJ \f'(e i0 ,t)r 



Now we turn to J\ 



2tt 2n 

J 1= [ Re (l + e ia ^] Re — [ — 4r 
J V f'J 27r7 |/'(e"», 



2e l V 



id Act 



Here we change the order of integration and get 



2tt 



2tt 







|/'(e*V)| 2 2tt 



Re — / Re 1 + e ^a ±- 



^)| 2 (e ie -e ia )' 



ia f"\ 2e»e* 



/' J (e ie - e 



ia\2 



da d6. 



Integrating by parts we obtain 

2tt 



2tt 



Jl 



1 



—I 



|/'(e^,t)| 2 2tt J da 



d 



Re — / —Re 1 + e iaJ — 



f"\ e ia + e i9 



f I e ia - e* 



:da dd. 
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The Cauchy formula gives 

2vr 



Jl=/ ^F Re [— + — -{— ] ]ie ' 



or 



j -/^ Re (K i+£ ?) 2+e2% -5i^ 

x 

where Sf is the Schwarzian derivative. 

Summing up all these integrals we come to the conclusion that 

d„„ I 1 „. /V, . e"f"\ , 2,»c . 1 



We observe that 

Re (l + ^) 



x(M) 



l/'(e*,*)l 

is the curvature of the boundary d£l(t) at the point f(e l6 ,t). Integration by parts implies 

' ^Im ( 1 + ^ 1 1 (W 



2n 



\f'(e i6 ,tW\ V /' 



So 



2 7 |/'(e^,t)| 2 V 2 V /' / 7 2 



2?r 2?r 

>, = |/^ (M)dS -3/_L_( Im ( 1 + ^ 

o o 

2vr 2vr 

+ /i?ra*(^' + S*-/i7<w" 

o o 

2tt 2tt 2 

= ( § + §)/^_§/^(m(i + *£)),» 


2vr 2tt 

+ / FW Re ( (1 + y° s < + <M>) * - / \nkw M - 
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Repeating this step we get 



2?r 2tt 



dO 



2vr 



7 l/'(e» 





n ■ 







|/'(e» f)l 2 



at n-th iteration. Taking limit as n — > oo, we finally obtain 

2tt 2tt 2tt 



| S [/] = 2/^(M)^ + / ^J-^Re^)^-/ 





as claimed in the statement of the theorem. □ 
In the simplest case of the circular evolution f((,t) = V%i( we have ^S[f] = J^£[/] = f ■ 

4. Parametric manifold Diff S 1 /Rot S 1 

As it has been mentioned in Section 2, starting with a smooth (C°°) initial boundary 
dQ(0) the classical evolution of Laplacian growth is given by domains Q(t) with smooth 
boundaries dQ(t) as long as the classical solution exists. Our aim now is to give an embed- 
ding of this evolution into the parametric Kirillov's space Diff S 1 /Kot S 1 . 

We denote the Lie group of C°° sense preserving diffeomorphisms of the unit circle 
S 1 = dU by Diff S 1 . Each element of Diff S* 1 is represented as z = e 1 ^ with a monotone 
increasing, C°° real- valued function (f)(9), such that (f>(9 + 2ir) = (f)(9) + 2ir. The Lie algebra 
for Diff S 1 is identified with the Lie algebra Vect S 1 of smooth (C°°) tangent vector fields 
to S* 1 with the Poisson - Lie bracket given by 

[01,02] = 0102 - 0201 • 

Fixing the trigonometric basis in Vect S 1 the commutator relations take the form 

r n „ n n — m.. n + m . 
[cos n9, cos mflj = — - — sin [n + m)U H — sin (n — mjfl, 

[sin n#, sin md\ = sin (n + m)9 H sin (n — m)9, 



n - 




m 




2 




rn 




n 




2 




rn 




n 


2 



[sin n9, cos m#] = "" cos (n + m)0 — n ~|~ m cos ( n _ 

There is no general theory of infinite dimensional Lie groups, example of which is under 
consideration. The interest to this particular case comes first of all from the string theory 
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where the Virasoro (vertex) algebra appears as the central extension of Vect S 1 (see Section 
7). The central extension of Diff S 1 is called the Virasoro-Bott group. Entire necessary 
background for the construction of the theory of unitary representations of Diff S 1 is found 
in the study of Kirillov's homogeneous Kahlerian manifold M = Diff S^/Rot S 1 , where 
Rot S 1 denotes the group of rotations of S 1 . The group Diff S 1 acts as a group of trans- 
lations on the manifold M with Rot S 1 as a stabilizer. The Kahlerian geometry of M has 
been described by Kirillov and Yuriev in JH]- The manifold M admits several represen- 
tations, in particular, in the space of smooth probability measures, symplectic realization 
in the space of quadratic differentials. Let A stand for the class of all analytic regular 
univalent functions / in U normalized by /(0) = 0, /'(0) = 1. We will use its analytic 
representation of M based on the class A of functions from A which being extended onto 
the closure U of U are supposed to be smooth on S 1 . The class A is dense in A in the local 
uniform topology of U. 

Let E stand for the class of all univalent regular maps in the exterior U* of the unit disk 
U normalized by g(() = ci( + cq + ci^ -1 + . . . which are smooth on S 1 . Then, for each 
/ G A there is an adjoint map g G E such that C\f(U) = g(U*). The superposition g^ 1 o / 
restricted to S 1 is in M. Reciprocally, for each element of M there exist such / and g. 
Observe that a piece- wise smooth closed Jordan curve is a quasicircle if and only if it has no 
cusps. So any function / from A has a quasiconformal extension to U*. By this realization 
the manifold M is naturally embedded into the universal Teichmuller space T. However, 
defined as a complex Banach manifold the Teichmuller space T requires additional Hilbert 
manifold structure to assure the embedding M — > T to inherit the Kahlerian structure of 
M. This has been done by Takhtajan and Teo in |SBl _ jlDl- The Kahlerian structure on M 
corresponds to the Kahlerian structure on T given by the analogue of the Weil-Petersson 
metric. 

The Goluzin-Schiffer variational formula lifts the actions from the Lie algebra Vect S 1 
onto A. Let / G A and let z/(e je ) be a C°° real-valued function in 9 G (0, 2ti] from Vect S 1 
making an infinitesimal action as 9 i— ► 6 + Tv(e %6 ). Let us consider a variation of / given 



Kirillov and Yuriev [15] . |2U] have established that the variations S v f(Q are closed with 
respect to the commutator and the induced Lie algebra is the same as Vect S 1 . Moreover, 
Kirillov's result J7j states that there is the exponential map Vect S* 1 — > Diff S 1 such that 
the subgroup Rot S 1 coincides with the stabilizer of the map f(() = ( from A. 



by 



(16) 
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5. Semigroups of conformal maps 



The basic ideas that we use in this section come from the development of Lowner's 
parametric method that emerges at a seminal Lowner's paper [23|. Lowner was first who 
proposed to use Lie semigroups of conformal maps to obtain an evolution equation for 
conformal maps. His ideas have been furthered by many authors among whom we mention 
Pommerenke (3111110] as a general reference, and Goryainov's works [TUt [TT] especially closed 
to our consideration, one also may see [S3 EZ] • 

We consider the semigroup Q of conformal univalent maps from U into itself with com- 
position as the semigroup operation. This makes Q a topological semigroup with respect 
to the topology of local uniform convergence on U. We impose the natural normalization 
for such conformal maps: $(£) = (3( + b 2 ( 2 + . . . , ( £ U, > 0. The unit of the semigroup 
is the identity. Let us construct on Q a one-parameter semi- flow $ T , that is, a continuous 
homomorphism from R + into Q, with the parameter r > 0. For any fixed r > the element 
$ T is from Q and is represented by a conformal map r) = /3(r)£ + b2(r)( 2 + . . . from 
U onto the domain $([/, r) C U. The element $ T satisfies the following properties: 

• $° = id; 

• $ r+s = $($(C, t), s), for r, s > 0; 

• r) —>■ £ locally uniformly in U* as r — > 0. 

In particular, /3(0) = 1. This semi-flow is generated by a vector field v(() if for each 
C £ U the function w = $(C, t), r > is a solution of an autonomous differential equation 
dw/dr = v(w) with the initial condition w\ T= o = C- The semi-flow can be extended to 
a symmetric interval (—t,t) by putting $~ T = $ -1 (£,r). Certainly, the latter function 
is defined on the set $(Z7, r). Admitting this restriction for negative r we define a one- 
parameter family $ T for r G (— £, t). 

For a semi-flow $ T on there is an infinitesimal generator at r = constructed by 
the following procedure. Any element $ T is represented by a conformal map r) that 
satisfies the Schwarz Lemma for the maps £7 — ► U, and hence, 



where the equality sign is attained only for $° = id ~ 0) = £. Therefore, the following 
limit exists (see, e.g., jTUHUllMl) 



Re 



c 



< 



$(C,t) 

c 



<i, Cef/, 



d<S>(C,T) 
dr 



T=0 



lim Re 

T^0 



Rc 



C 



< 



and the representation 



dr 



T=0 



(p(0 
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holds, where p(() = p + pi( + . . . is an analytic function in U with positive real part, and 

d(3(r) 



= -Po- 

r=0 



^ 3r 
In ^2] it was shown that $ T is even C7°° with respect to r. The function — Cp(C) is an 
infinitesimal generator for $ T at r = 0, and the following variational formula holds 

(18) $(C,t) = C-tCp(0 + °(t), ^(r) = l-rp + o(r). 
The convergence is thought of as local uniform. We rewrite (|18|) as 

(19) $((, r) = (1 - rp )( + r C(-p(C) + Po) + o{r) = (3(r)( + r ((-p(() + p ) + o(r). 

Now let us proceed with the semigroup Q G Q of elements from Q represented by univalent 
maps $ smooth on S . By the variation of the identity in A given by the formula (|16|) we 

get 

$(C,r) C 2 [_dH_, M f x 



/3(r) 27ri y w(w — () 

s 1 

for v(w) from Vect S* 1 . Comparing with (|19|) we come to the conclusion about p(C) : 
(20) p(2)=po + _£_| 



iu(to - C) 



The constants po and the function v(w) must be such that Re p(z) > for all z £ U. 
We summarize these observations in the following theorem. 

Theorem 2. Let $ r be a semi-flow in Q. Then it is generated by the vector field v(() 

-Cp(C), 



-dw, 



w(w - C) 
s 1 

where v{e %e ) G Vect S* 1 , and i/ie holomorphic function p(Q has positive real part in U . 
This theorem implies that at any point r > we have 

??&li = _$(C,T)p(*(C,r)). 

6. Evolution families and evolution equations 

A subset $ i,s of < s < t is called an evolution family in if 

• = id; 

• $*> s = $*' r o $ r < s , for < s < r < i; 

• — > id locally uniformly in U* as t, s — > r. 
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In particular, if $ T is a one-parameter semi-flow, then $*~ s is an evolution family. We 
consider a subordination chain of mappings f((,t), ( G U, t G [0,to), where the function 
= a(t)z + ct2(t)( 2 + . . . is a analytic univalent map U — > C for each fixed t and 
/([/, s) C f(U,t) for s < t. Let us assume that this subordination chain exists for t in an 
interval [0, T). 

Let us pass to the semigroup Q. So now has a smooth extension to S 1 . Moreover, 
$' ,s — > id locally uniformly in C as t, s — > r. 

We construct the superposition / _1 (/(C, s )j*) f° r £ £ [0, T") , s < t. Putting s = t — r we 
denote this mapping by r). 

Now we suppose the following conditions for f((,t). 

(i) The maps f{(,t) form a subordination chain in [/, £ G [0,T). 

(ii) The map f(C,t) is holomorphic in U, f((,t) = a(t)( + a 2 (t)( 2 + where 
a(t) > and differentiable with respect to t. 

(iii) The map f((,t) admits a smooth continuation onto S 1 . 

The function <!>(£,£, r) is embedded into an evolution family in Q. It is differentiable 
with regard to r and t in [0, T), and &((,t, 0) = Moreover, £ = lim T _^ $(C> ^ r ) locally 
uniformly in [/ and t, r) is embedded now into an evolution family in Q qc . The identity 
map is embedded into a semi-flow $ r C ^ (which is smooth) as the initial point with the 
same velocity vector 

— ^ =-CKC.*), Ce^, 

ar r=o 

that leads to the Lowner-Kufarev equation 

(21) f = Cfp{C,t), 

(the semi-flow $ T is tangent to the evolution family at the origin). Actually, the differ- 
entiable trajectory f((,t) generates a pencil of tangent smooth semi-flows with starting 
tangent vectors —(p{(,t) (that may be only measurable with respect to t). 

Therefore, the conclusion is that the function f{(,t) satisfies the equation (|2ip where 
the function p((,t) is given by 

, . v , . C f u(w,t) , 
P(C, t) = Po {t) + t^t / / ' dw, 

s 1 

and has positive real part. The existence of po(t) comes from the existence of the subor- 
dination chain. One may assign the normalization to f(C,t) controlling the change of the 
conformal radius of the subordination chain by, e.g., e*. In this case, changing variables we 
obtain p = 1. 

Summarizing the conclusions about the function p((,t) we come to the following result. 

Theorem 3. Let f{(,t) be a subordination chain of maps in U that exists for t G [0,T) 
and satisfies the conditions (i-iii). Then, there are a real valued function p (t) > and a 
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real valued function t) G Vect S 1 , such that Re p((, t) > for ( G U , 

P(C, t) = Po(t) + ^ I U } W,t \ dw, C G 17, 



27ri 7 w(w — 

an<i f{(,t) satisfies the L owner- Kufarev differential equation \21\) in t G [0, T). 

Comparing (jlHj) and (j!4j) with this theorem we come to the following corollary. 

Corollary 1. Let f((,t) = a(t)( + c^WC 2 + • • • be a subordination chain of maps in U 
that parameterizes the classical Laplacian Growth that exists fort G [0,T). Then, under 
the notations of the previous theorem we have 

2vr 

a 1 f 1 7 / id \ 2 

Po = - = 7T / 77777F~rTT?" 6/ ' and u \ e >*) 



a 2W |f(e^,t)| 2 ' v ' y |/'(e* t)P" 

o 

In the view of this corollary we come to an interpretation of the variation of the loga- 
rithmic action as 

2tt 2tt 

±(S[f] + S[f}) = 2 J x 2 (6,t)d6 + Re J e™v(e> e ,t)S f d6, 
o o 

where v is a vector from the Lie algebra Vect S 1 tangential to Diff S' 1 /Rot S 1 at the unity. 

7. Connections with the Virasoro algebra 

In two dimensional conformal field theories |H] , the algebra of energy momentum tensor 
is deformed by a central extension due to the conformal anomaly and becomes the Virasoro 
algebra. The Virasoro algebra is spanned by elements = ( 1+k d, fceZ and c with efc+e_fe, 
where c is a real number, called the central charge, and the Lie brackets are defined by 

[e m , e n ] V ir = (n - m)e m+n + —m{m 2 - l)S n _ m , [c, L k ] = 0. 
The Virasoro algebra (Vir) can be realized as a central extension of Vect S 1 by defining 

[<pd + ca, ipd + cb] Vir = (0-0' - (j>'ip)d + y^O, 

(whereas [4>,ip] — fyty' ~ </>V)> where the bilinear antisymmetric form u(<f),ip) on Vect S 1 is 
given by 

2tt 


and a, b are numbers. This form defines the Gelfand-Fuks cocycle on Vect S 1 and satisfies 
the Jacobi identity. The factor of 1/12 is merely a matter of convention. The manifold M 
being considered as a realization A admits affine coordinates {c2, C3, . . . }, where Ck is the 
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k-th coefficient of a univalent functions / G A. Due to de Branges' theorem [Sj, M is a 
bounded open subset of {\c k \ < k + e}. 

Taking v = —ie lk0 , k > 0, we obtain the expressions for L k = 5 u f, f G A ~ M (see 
formula (JUJ)), as 

A> = C/'(C)-/(C), L k = c 1+k f. 

The computation of for k < is more difficult because poles of the integrant. For 
example, 

L_ x = /' - 1 - 2c 2 /, L_ 2 = ^ - 1 - 3c 2 + (c| - 4c 3 )/, 
(see, e.g., |TS] ) - In terms of the coordinates {c 2 , c 3 , . . . } on M 

oo oo 

L k = d k + + l)c n d k+n , L = ^2 nc nd n , 

n=l n=l 

for k > 0, where 9^ = <9/<9c fc+1 . 

Neretin |2E1 introduced the sequence of polynomials P k , in the coordinates {c 2 , C3, . . . } 
on M by the following recurrent relations 

L m {P n ) = (n + m)P n _ m + ^m(m 2 - l)5 n , m , P = Pi = 0, P fe (0) = 0, 

where the central charge c is fixed. This gives, for example, P 2 = | (03 — c 2 ), P3 = 2c(c4 — 
2c 2 C3 + c|). In general, the polynomials P k are homogeneous with respect to rotations 
of the function /. It is worthy to mention that estimates of the absolute value of these 
polynomials has been a subject of investigations in the theory of univalent functions for a 
long time, e.g., for |P 2 | we have \c$ — c|| < 1 (Bieberbach 1916 0j), for estimates of IP3I 
see |13| |2HJ El H2] • For the Neretin polynomials one can construct the generatrix function 

k=l 

where Sf(() is the Schwarzian derivative of /, Let v G Vect S 1 and v 9 be the associated 
right-invariant tangent vector field defined at g G Diff S 1 . For the basis v k = —ie lhe d, one 
constructs the corresponding associated right-invariant basis v 9 k . By {ip~ k } we denote the 
dual basis of 1-forms such that the value of each form on the vector v 9 k is given as 

i^k, K) = S k+n,0- 

Let us construct the 1-form Q on Diff S 1 by 

00 

Q = ^(Pfc o n)ip k , 
k=i 
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where it means the natural projection Diff S 1 — > M. This form appeared in E] in the 
context of the construction of a unitarizing probability measure for the Neretin representa- 
tion of M. It is invariant under the left action of S 1 . If / £ A represents g and v e Vect S 1 , 
then the value of the form Q on the vector v is 



o 

see [21 E] • So the variation of the logarithmic action given in Theorem 1 becomes 

2tt 





In this formula, we take into account the first coefficient, the conformal radius of the 
Laplacian evolution, that does not change the form Q. 



(i) One may conjecture that the formula ()22|) remains true for general smooth subor- 
dination evolution. 

(ii) The logarithmic action is clearly related to the universal Liouville action suggested 
by Takhtajan and Teo in |38lEMHUj . There must be possible to obtain the variation 
given in Theorem 1 by means of the variation of the universal Liouville action 
obtained in ■ This would be yet more interesting because the universal Liouville 
action is defined for contours without any smoothness hypothesis. 

(iii) Another interesting question is whether it is possible to make regularization of 
the proper Liouville action integral based on the Poincare metric by the boundary 
distortion by a univalent function. 

(iv) A deeper task is concerned with the Laplacian growth and its embedding into the 
Whitham-Toda hierarchy. The extended Toda hierarchy (see jH]) admits a non- 
abelian algebra of infinitesimal symmetries isomorphic to half of the Virasoro alge- 
bra [7j. It would be interesting to reveal the connections between the Hamiltonian 
approach through the Toda hierarchies and the action approach suggested in the 
present paper. 

(v) The multiply connected Laplacian growth is a natural way of generalization of all 
these results. 
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